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Supplementary Experimental Methods 

Strains construction 

We engineered a two-color fluorescent reporter strain {Escherichia coli) to measure 
the activity of the marRAB promoter, following the work by Miyashiro & Goulian 
(2007) \ The strain (named as IE01) contains a chromosomal copy of the yfp gene 
(Yellow Fluorescent Protein) under the control of the marRAB promoter, and the cfp 
gene (Cyan Fluorescent Protein) expressed with a constitutive tetA promoter at the 
attachment sites in E. coli of lambda and HK022 phages, respectively. 

To construct the YFP reporter, we amplified the promoter region of the 
marRAB operon (1616883-1617144) by PCR. The primers used were: 5'- 
TAGC AGAATTCCGGCAGCAAC (forward), and 5 ' - 

GCATAGGGATCCTGGCAAGTAAT (reverse). The PCR product was then cloned 
it into a plasmid (pTM74; Miyashiro & Goulian, 2007) that contains a multicloning 
site upstream of the promoterless yfp gene to yield the plasmid pJFPOl (this work). 
This plasmid was integrated as single copy into the lambda phage attachment site of 
E. coli MG1655 using the helper plasmid plnt-TS (Weiss et aL, 1999). The residual 
cat gene linked to the marRAB promoter was replaced by kan by electroporating a 
PCR product created from the template plasmid pKD4 into the appropriate strain 
carrying plasmid pKD46 to yield the strain JFP01 (this work). The two-color reporter 
strain JFP02 was created by moving the previous marker into strain TIM64 (a strain 
derived from a MG1655 that constitutively expresses CFP; TIM64 does not present 
any antibiotic resistance marker; Miyashiro & Goulian, 2007) by PI transduction. 
Removal of the kan marker using plasmid pCP20 originated E. coli strain IE01. This 

1 We thank Tim Miyashiro for strains, and Mark Goulian for strains and experimental advice. 
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was verified by PCR. In addition, E. coli strain IE02 was constructed by deletion of 
the rob gene in IE01 with the application of a Datsenko & Wanner (2000) knockout 
protocol followed by the removal of the kan marker using plasmid pCP20 (Datsenko 
& Wanner, 2000). The primers used were: 5'- 
ATATCCCAATGGCATCGTAAAGAACATTTTGAGGCATTTCAGTCAGTTGC 
GCTGGAGCTGCTTCGAA (forward), and 5'- 

ATGAACCTGAATCGCCAGCGGCATCAGCACCTTGTCGCCTTGCGTATAAT 
ATGAATATCCTCCTTAG (reverse). The sequence of the resulting strain was 
verified by PCR. 

We also constructed a strain deleting the marB gene (named as DB01) 2 . This 
was done from strain IE01 using a variant of the Datsenko & Wanner (2000) 
knockout protocol, with the insert {kan marker with marB flanking regions) obtained 
by PCR from a KEIO collection strain (Baba et al. 9 2006). Strain IE01 containing the 
pKD46 plasmid was electroporated with this insert in presence of arabinose, allowing 
recombination in the marB locus. The primers used to obtain the insert were: 5'- 
GACCAATATGC AGGGCGAATCG (forward), and 5 ' - 

ATGTATTTGGCTTGCGGTGG (reverse). Correct substitution of the marB gene by 
kan marker was checked by PCR and sequencing. 

Culture media and reagents 

Medium LB was always used for overnight cultures. Minimal medium M9 (M9 salts 
lx, MgSCU 2 mM, CaCb 0.1 mM, glucose 0.4%, casamino acids 0.05%, vitamine Bl 
0.05%)) was used to grow cells during characterization experiments. To induce the 
mar circuit, we used different concentrations of Salicylate (Sigma), resulting in a 

2 

We thank Jeronimo Rodriguez-Beltran for practical assistance. 
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gradient ranging from 0 to 5 mM, and also different concentrations of cAMP (Sigma), 
resulting in a gradient ranging from 0 to 20 mM. When appropriate, we used 
kanamycin at a concentration of 50 |ig/mL. 

Quantification of fluorescence in a cell population 

Experiments of induction were carried out to study the dynamics of the system. 
Cultures (volume of 2 mL) inoculated from single colonies (three replicates) were 
grown overnight in LB medium supplemented with glucose 0.4% at temperature of 37 
°C and shaking of 170 rpm. Cultures were then diluted 1:200 in M9 minimal medium 
(10 |iL of culture into a final volume of 2 mL) and were grown for 2 h at temperature 
of 37 °C and shaking of 170 rpm. Cultures were then used to load the multiwell plate 
(Thermo) with final volumes of 200 jiL. Per well, we added directly a volume of the 
culture, a volume of salicylate stock, and a volume of cAMP stock (when 
appropriate) 3 . 

The set up for the Victor X2 was as follows. OD600 measured with 
absorbance filter of 600 nm (0.5 s for reading). YFP measured with excitation filter of 
497/16 nm, and emission filter of 535/40 nm. CFP measured with excitation filter of 
434/17 nm, and emission filter of 479/40 nm. For both YFP and CFP, we took 8000 
for energy lamp, and 0.1 s for counting. Temperature was at 37 °C. The program 



We also used mineral oil for characterization. Cultures were grown as previously indicated, and then were used 
to load the multiwell plate (Thermo) with final volumes of 250 ^L. Per well, we added directly 200 of culture 
(counting the volume of inducers), and 50 of mineral oil (Sigma). However, we found that the mar circuit does 
not respond to salicylate in anaerobic conditions. 



4 



started first with OD600, then YFP, and finally CFP, followed by 30 s of shaking in 
orbital mode. Then it waited for 5 min, and it started again. 

Analysis of fluorescence data 

We collected time-course data of fluorescence (YFP and CFP) and absorbance. 
Background values of absorbance and fluorescence, which corresponded to M9 
minimal medium, were subtracted to correct the signals. The normalized fluorescence 
(for both YFP and CFP) was calculated as the ratio of fluorescence and absorbance. 
Similar values of normalized fluorescence were reported for MG1655 cells and for 
M9 minimal medium, which indicated that the auto-fluorescence of cells was 
negligible in this case. The growth rate of cells was calculated as the slope of the 
linear regression between the log of corrected absorbance and time in exponential 
phase. Time-dependent promoter activity, defined as the instantaneous production rate 
of normalized YFP fluorescence (magnitude per cell), was calculated for each time 
point using the derivative of the normalized fluorescence. Promoter activity in steady 
state, defined as the stationary production rate of normalized YFP fluorescence 
(magnitude per cell) in exponential phase, was calculated as the average over time 
(for t > 2 h) of normalized fluorescence times growth rate. The error associated to that 
measure was obtained by calculating the standard deviation over replicates in all time 
points, then squaring all these deviations to average them, and finally getting the root 
square. Data were analyzed with Matlab (Math Works). 

Quantification of mRNA levels by qRT-PCR 

Cultures (volume of 2 mL) inoculated from single colonies (three replicates) were 
grown overnight in LB medium at temperature of 37 °C and shaking of 170 rpm. 
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Cultures were then diluted 1:200 in M9 minimal medium (10 jiL of culture into a 
final volume of 2 mL) and were grown for 2 h at temperature of 37 °C and shaking of 
170 rpm. Then, aliquots of 1.5 mL were quickly pelleted by centrifuging 2 min at 
13,000 rpm and resuspended in 50 ^iL of TE buffer (10 mM Tris-HCl, pH 8.0, 1 mM 
EDTA). Bacteria were broken by adding 50 \xL of a mix 1:1 phenol: chloroform (pH 
8.0) and vortexing thoroughly. Bacterial RNA from each sample was recovered in the 
aqueous phase by centrifuging 5 min at 13,000 rpm. A subsequent re-extraction was 
done by adding 50 |iL chloroform followed by centrifugation of 5 min at 13,000 rpm, 
recovering the aqueous phase. Samples were then passed through a silica-based, 
DNA-clean column (Zymo) and were eluted in 8 jiL of TE buffer. DNase I 
(Fermentas) was added (1 U) for 1 h at 37 °C, which was then thermally inactivated 
(30 min at 65 °C) by adding EDTA. Total RNA eluted was quantified in a NanoDrop. 

One-step SYBR PrimerScript RT-PCR Kit II (Takara) was used for 
detection, following the Kit protocol for preparing the reaction volumes. 16S rRNA 
was used as housekeeping gene to normalize RNA quantity in each reaction. The 
primers for qRT-PCR were synthesized by IDT for mar A and Sigma for 16S rRNA. 
The primer sequences for amplifying mar A were taken from (Viveiros et al. 9 2007) 
and for 16S rRNA from (Carothers et al. 9 2011). Reactions in triplicate were carried 
out using a Step One Plus Real-Time PCR System (Applied Biosystems). The thermal 
cycling program for amplification was 5 min at 42°C, 10 s at 95°C, and 40 cycles of 5 
s at 95°C and 34 s at 60°C (Shuttle PCR); following with default melting curves 4 . 



We thank Eszter Majer for help with the RNA extraction and Anouk Willemsen for help with the qRT-PCR. 
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Quantification of fluorescence in single cells 

Experiments of induction were carried out to study the heterogeneity of the dynamic 
response. Culture (volume of 2 mL) inoculated from a single colony was grown 
overnight in LB medium supplemented with glucose 0.4% at temperature of 37 °C 
and shaking of 170 rpm. Culture was then diluted 1:200 in M9 minimal medium (10 
|iL of culture into a final volume of 2 mL) and was grown for 4 h at temperature of 37 
°C and shaking of 170 rpm. The different cultures diluted 1:10 were then used to load 
the agarose pads. 

Agarose pads were prepared with a volume of 5 mL of M9 minimal medium 
and 0.075 g of agarose. It was dissolved by vortexing and microwaving. The pads 
were then allowed to solidify for about 1 h at room temperature before seeding 
bacteria. 2 |iL of each culture were then used to load the pads. They were kept for 
about 15 min at room temperature so that cells can be absorbed into the agarose. Just 
before characterization with the microscope, 2.5 \iL salicylate from a solution 0.1 M 
was used to induce cells in solid medium, having estimated the volume of the agarose 
pad in 50 jiL (resulting concentration of salicylate about 5 mM). 

In each pad, fields with an adequate initial density of cells were chosen, the 
first photo was taken, and salicylate was added. Photos were taken for each field 
every 12 minutes. 

Agarose pads were monitored in an inverted microscope Axiovert200 (Zeiss) 
with objective 100X/1.45 oil Plan-Fluar at temperature of 30 °C. The microscope was 
equipped with a digital camera C9 100-02 (Hamamatsu), a Xenon lamp XBO 75W/2, 
an optical filter changer Lambda 10-2 (Sutter), a motorized stage (Marzhauser), and a 
temperature controller (Cell Observer-Zeiss). Moreover, the microscope was 
automated by the commercial software MetaMorph (Universal Imaging). Cell images 
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were acquired from the bright-field and fluorescence channels. We used the 
fluorescence filters yellow FP (490-510 nm, 510-560 nm) and cyan FP (426-446 nm, 
460-500 nm). 

Analysis of single cell images 

Microscopy photos were segmented and analyzed using the EBImage package for R 
from Bioconductor. Segmentation allowed us to identify sets of pixels belonging to 
individual cells, and measure the average apparent YFP and CFP intensities for each 
of these sets. The ratio YFP/CFP was used as proxy for system activity. 

Quantification of fluorescence by flow cytometry 

Cultures (volume of 2 mL) inoculated from single colonies were grown overnight in 
LB medium without and with 5 mM salicylate at temperature of 37 °C and shaking of 
170 rpm (to have uninduced and pre-induced cultures). These cultures were 
subsequently diluted 1:200 in 2 mL of M9 minimal medium with different 
concentrations of salicylate (0, 0.05, 0.25, 0.5, 2.5 and 5 mM). After incubation for 4h 
at the same conditions, cells were spun down and resuspended in PBS buffer. YFP 
(525 nm) intensity was measured in the flow cytometer Gallios (Beckman-Coulter). 
YFP distributions were obtained with a standard flow cytometry analysis (software 
Kaluza). 



8 



Supplementary Mathematical Methods 

Empirical model of the dynamic response 

Our reporter protein (YFP) monitors the promoter activity of the marRAB operon 
(Pmar)- Then, we can write 

^-[YFP] = Il mar -ii[YFP] (SI) 
at 

where is the cell growth rate and H mar the activity of promoter P mar . Empirically, we 
proposed the following exponential model 

[YFP} m ^ 1 

where A and m are two parameters that control the dynamics of the response, and 
[YFP], is the concentration in steady state. We fitted with our own experimental data 
by nonlinear regression the parameters A and m. Combining Eqs. (SI) and (S2), 
promoter activity was calculated as 

The time to reach the 50% of the steady state value (tso) upon induction with 
salicylate reads 



-1, 

' 50 = T 



i- 1 



2 l < 



(S4) 



1 d[YFP] 

We also considered the response speed, here defined as and evaluated 

[YFP] dt 

at tso. This gives ln(2) A for Eq. (S2) and m » 1 . In case of a non-regulated gene, 

whose dynamics is governed by (l-e'^j, the response speed is directly \x. In 

addition, differentiating Eq. (S3), we calculated the time where promoter activity is 
maximal (r), given by 
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1, 

r = —In 
A 



(S5) 



We used this model to fit either population or single cell experimental data. 



Bottom-up mathematical model 

We constructed a system of ordinary differential equations (ODEs) by knowing the 
topology of the circuit. MarR, MarA and MarB form an operon controlled by 
promoter P mar (Alekshun & Levy, 1997; Chubiz et al. 9 2012). MarR represses P mar , 
which can be modulated by salicylate (Cohen et aL, 1993). MarB also represses P mar 
(Vinue et al. 9 2013), whereas MarA activates it (Martin et ai, 1996). On the other 
hand, Rob is controlled by promoter P ro b, and it activates P mar (Alekshun & Levy, 
1997). MarA and Rob repress P ro b (Schneiders & Levy, 2006). Finally, in our system 
YFP models a downstream gene controlled by promoter P mar . Although MarR, MarA 
and MarB are transcribed from the same promoter, the corresponding protein 
expressions may be different each other due to distinct translation rates. By analyzing 
the 5' untranslated regions of MarR, MarA and MarB with RBS calculator (Salis et 
aL, 2009), considering the 30 nucleotides upstream and the 7 nucleotides downstream 
of the start codon, we found that translation rates of MarA and MarB are about 30- 
fold and 20-fold, respectively, higher than the translation rate of MarR. This is in tune 
with previous experimental observations (Martin & Rosner, 2004). In addition, 
promoter P mar is regulated by CRP-cAMP. Therefore, we could write 
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^[MarA] = /M^ - + 6)[MarA] 
at 

^[MarR] = U nmr - f j[MarR] 
at 

^[MarB] = /3*U mar - ^[MarB] (S7) 
dt 

^[Rob] = U rob -fj[Rob] 
dt 

^-[YFP] = U mar -^YFP], 
at 

where is the cell growth rate, 8 the degradation rate of MarA {8 » noting that 

MarA is quickly degraded by protease Lon (Griffith et al. 9 2004), (5 (and /f ) the fold 

increase of MarA (MarB) translation rate, and H mar and Il ro b the activity of promoters 

P mar and Prob, respectively. 

The equations for H mar and H ro b, knowing that MarA and Rob act as 

monomers whereas MarR as a dimer (Martin et ai, 1996), could be approached by 

Hill functions (Bintu et ai, 2005). However, we knew that the effect of MarA on P ro b 

is not observable in physiological conditions (Chubiz et aL, 2012), so we neglected 

this regulation to simplify the system of equations. We also demonstrated in this work 

by qRT-PCR that MarB does not repress P mar in presence of salicylate, whereas MarB 

indeed represses P mar in absence of salicylate. This also led us to eliminate MarB 

from the model. Therefore, it turned out 

n =n 1 + f a [MarA]/K A + f b Q([Rob],[cAMP]) 1 
mar 0 l + [MarA]/K A+ Q([Rob],[cAMP]) 1 + {[MarR fr J/K R f ' 

where Ka 9 and K R are the effective dissociation constants for transcription regulation, 
and f a and the activation fold changes. IIo is the basal protein synthesis rate, and Q a 
regulatory function that accounts for the activation of Rob and CRP-cAMP. 
Phenomenologically, and assuming competitive binding between the three 
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transcriptional activators (MarA, Rob and CRP-cAMP) as their operators overlap, we 
can write 

Q([RobUcAMP]) = ([Rob]/K B ) + ([cAMP]/K c ), (S9) 
where Kb and Kc are the effective dissociation constants for transcription regulation. 
This model can explain that, in presence of Rob or cAMP, MarA does not increase 
significantly the occupancy of polymerase at the marRAB promoter, which agrees 
with previous experimental and theoretical data (Wall et al, 2009). In addition, we 
have 

* l + ([Sal]/e s ) Vs 

where Os is the effective dissociation constant between salicylate (Sal) and MaR, Vs 
the Hill coefficient, and a the minimal fraction of free MarR. 



Note on the negative autogenous control of the system 

In addition to MarR and MarA, the marRAB operon contains another open reading 
frame coding for a small periplasmic protein. Recently, MarB has been identified as a 
repressor of the marRAB promoter (Vinue et ai, 2013), establishing a new putative 
negative feedback into the system. To study this novel regulation on the dynamic 
response, we knocked out the marB gene. A quantification of transcripts revealed that 
MarB is not able to repress the expression of the marRAB operon in presence of 
salicylate, although MarB does repress the operon in absence of it (Fig. S2N), which 
is tune with previous results (Vinue et al, 2013). We confirmed by fluorescence 
assays the observation that MarB does not repress the marRAB promoter upon 
induction with salicylate, even for moderate concentrations of this inducer (Fig. 
S2M). While MarB does not impact the dynamic range, we observed that the kmarB 
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system exhibits a slower response to salicylate. Our results suggested that the 
bacterial cell might exploit MarB as an enhancer of MarR to tightly repress the 
operon. They also revealed that Rob and MarB are regulatory elements that contribute 
to a fast mounting of the mar phenotype. 



Simplification of the mathematical model 

Our bottom-up mathematical model can be simplified for a better analysis of the 
dynamic response. By noting x = [MarA]/K R , y = [MarR]/K R , y 0 = [MarRf ree ]/K R (the 
concentration of MarR is equal to the one of YFP), p =f a =fb (assuming for simplicity 
equal induction by the three transcriptional activators), Jto = Uq/Kr, and also k = 
KrIKa, we could write a simplified system of ODEs. We also denoted by z ([Rob]/K B ) 
and c ([cAMP]/Kc) the normalized concentrations of Rob and cAMP, respectively 
(and we took K B = K A ). Thus, we obtained 

— = I3fc mnr - 6x 
at 

^--x^-W (Sll) 
at 



1 + p(kx + KZ + c) 1 



L mar JL () , t 2 ' 

l+KX + KZ + C 1 + Jo 

where MarA could be approached to a quasi-steady state (x^Jt mar , a function of 
time). And also we had 

y 0 l + a{[Sal]/d s y s 

y l + ([Sal]/d s ) Vs 

In case of maximal induction of the system with salicylate, a modulates the 
regulatory role of MarR. To obtain dimensional parameters, see Table S 1 for values 
of Kr and Ka. 
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Analytical solutions of the model I: temporal behavior 

In this section, we considered that a very strong the repressor acted on the system, 
MarR (i.e., y » 1). We also assumed that the system was induced with high levels of 
salicylate, so that yo = ay, and that the activation term was simply reduced to the fold 
change (p). This allowed us to simplify the model to just one ODE, given by 

dy Jt 0 p 



2 2 



dt ay 



-fjy. (S13) 



/ \ 1/3 

The steady state of this ODE was = — ^-y . This has the following solution 

\fia ) 

(after integration by separation of variables) 

A = (i_ e -W /3 . ( si4) 
y 

J 00 

However, this analytical solution gave an over-estimation of the dynamics when 
solving numerically the ODE. 

We then considered a situation at short times (i.e., ay « 1) to write 

The steady state of this ODE was y 00 =—. This has the following solution (after 

a 

integration by separation of variables and partial fraction decomposition) 
— = tanh( jt 0 pat) . (S16) 

3^00 

We verified with numerical simulations that this analytical solution is a good 

approximation for t <t 50 . With this model, we calculated the response time of the 

ln(3) .. , . 1 . 

system as r 50 = (i.e., y(r 50 ) = -yJ. 

2ji 0 pa 2 
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Analytical solutions of the model II: dynamic range 

Our model to describe the dynamics of a self-repressed, self-activated operon (y), 
implemented with two regulatory genes (repressor and activator), can be written as 

4_ * 0 i+P^'y (S17) 

dt l + (a s y) 2 l + K'y 
where as describes the effect of the inducer molecule (Eq. SI 2), which acts on the 
repressor (1 > a s > a). Here, MarA was assumed to be a stable protein (6 = fj). This 
assumption leads to x = /Jy, and does not change the dynamics of our protein of 
interest, MarR (y). This model can be solved in steady state for a strong repressor (y 
» 1) and high activation fold (p » 1) as 



( \ 1/3 

2 



*a s - 2l \ (SI 8) 



For saturating levels of salicylate, we estimated with this model an output dynamic 
range of R out = <x 2/3 (i.e., the ratio between the highest and lowest values of when 

varying salicylate), and an input dynamic range of R[ n = 9 2lnin (i.e., the ratio between 
the salicylate values at which we have y^=mm(y QO ) + 0.\lSy QO and 
y^ =max(j oo )-0.1Aj oo , where Ay w = max(y J - min(y J ; n in is an effective Hill 
coefficient; see also Goldbeter & Koshland, 1981). For a = 0.05, it turned out R out = 
7.37 and R in = 28.32 ( n in = 1 .3 1). 

For a constitutively expressed regulator (y), the dynamics of its regulated 
operon ( y ) can be written as 

f-rPV** (S19) 

dt l + (a s y) 

where y is constant. The steady state solution for y is straightforward to obtain. 
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Moreover, a model to describe the dynamics of a self-repressed gene (y), with 
equal production rate as before, can be written as 

f=-fV^- (S20) 

dt l + (a s y) 

This model can be solved in steady state to obtain the same expression as before. 

Finally, a model to describe the dynamics of a self-repressed gene (y) that 
becomes self-activator in presence of the inducer can be written as 

dy l+p((l-a s )yf 

— = ^o — ; 7^-7-, '-—2-w, (S21) 

dt i + ( asV ) 2 + ((i_ as ) v ) 

where we have assumed competitive binding between the repressor and activator. 
This model can be solved in steady state, also for a strong repressor (y » 1) and high 
activation fold (p » 1), as 

y„«^d{\-a s ), (S22) 
\x 

where 6 is the Heaviside function. 



Signal transduction and de-repression 

Our model also explains why the translation rate of MarR is limited. For the repressor 
to work as an efficient sensor, it has to reach an appropriate concentration. First, the 
concentration has to be higher than the effective dissociation constant (Kr, interaction 
protein-DNA). In this way, the regulator can exert its repressive action on the 
promoter (in absence of inducer). In the case of MarR, K R is of the order of nM (Table 
SI), which allows having [MarR] > K R (lower bound). Second, at maximal 
concentration of inducer (Sat), the promoter has to be de-repressed. This imposes the 
condition [MarRf ree ] < K R . In the marRAB operon, the environmental molecule 
(antibiotic or salicylate) does not interact directly with the repressor (contrary, for 
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example, to the lac operon). Instead, Cu is responsible for MarR oxidation (Hao et 
al 9 2014), and the intracellular concentration of this cation is up to 10 |iM (Rensing & 
Grass, 2003). This clearly imposes a maximal concentration of repressor. Knowing 
the oxidation curve of MarR by Cu 2+ (Fig. S1G, and taking the effective oxidation 
constant as K ox « 1), it turns out [MarR] < (K R [Cu 2+ ] 2 ) 1/3 (upper bound). Therefore, 
we propose that an inefficient translation rate of MarR might have evolved to 
appropriately deal with copper signaling. 



Apparent transfer function and stress cross-talk 

T-rmin rrmax/r SdH / K Y S 

The Hill-like model X\ mar = ^ — J — ^— was adjusted independently to 

l+([Sat]/K s ) ns 

the experimental data of induction of the wild-type and Arob systems, obtaining the 
corresponding apparent transfer functions ([Sal] is the concentration of salicylate in 
mM). We got = 2,260 (+/- 746) AU, n™* = 20,000 AU (fixed), K s = 0.88 (+/- 
0.12) mM, and n s = 1.28 (+/- 0.19) for the wild-type system, and = 1,089 (+/- 
1011) AU, ITZ = 12,500 AU (fixed), K s = 0.72 (+/- 0.19) mM, and n s = 1.45 (+/- 
0.48) for the Arob system. With the inferred sigmoidal models, we calculated the 
input and output dynamic ranges: R in = 20- 30 and R out = 9 - 1 1 for both systems. 

Apart from salicylate, the marRAB operon can also be induced by CRP- 
cAMP. We predicted computationally and verified experimentally the two- 
dimensional transfer function of the system. Then, we defined a degree of stress-cross 



talk as 



\cAMP 



\YFP\ 



-1 



\cAMP = 0 
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Stochastic modeling 

Using a Langevin approach, our ODE-based mathematical model can be extended to 
account for the inherent stochasticity of biological systems. Because in bacteria noise 
has two components (intrinsic and extrinsic) and is predominantly generated at the 
transcription level (Swain et al. 9 2002), we considered a stochastic process (where 

jji is the inverse of the correlation time) with statistics ^(f)) = 0 and 

^(f')^(^+0) = ^e~^. Thus, and having assumed similar mRNA degradation and 
MarR translation rates, we can write 



where q is the extrinsic noise magnitude. This system can be solved numerically by 
following the method described by Rodrigo et al (201 1). 

Note on the potential excitability of the system 

The regulatory architecture of the marRAB operon is suitable to present excitable 
behavior (even to produce sustained oscillations) in the non-induced state, due to the 
interlinked, delay-based positive and negative feedbacks (Hasty et ai, 2002; Guantes 
& Poyatos, 2006). For this architecture, oscillations appear when a clear separation of 
time scales between activator and repressor dynamics exists. This implies two 
conditions: /) the activator degradation should be stronger than the repressor one (in 
our case, MarA is quickly degraded by a protease, whereas MarR is diluted due to cell 
growth rate), and ii) the activator translation rate should be also stronger than the 
repressor one (here, the RBS of MarA is stronger than the one of MarR). However, in 




(S23) 
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the mar system, the activator (MarA) works as a monomer, which prevents this 
behavior even having much faster dynamics (see in Fig. S5B how oscillations arise 
when high cooperativity is introduced for both regulators). Moreover, when 
competitive binding between MarA and MarR (in particular, that MarA prevents 
MarR binding; Garcia-Bernardo & Dunlop, 2013) is considered, stochastic pulses can 
appear, even with a monomeric activator (Fig. S5E). However, previous experiments 
support relatively independent action of MarA and MarR (Martin & Rosner, 1997), 
aside from their operators do not overlap. Furthermore, in our single cell experiments 
without salicylate, we did not observe pulses in YFP, neither with our predictive 
model. 
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Supplementary Figures and Figure Legends 




time (h) time (h) 



C D 




Figure SI. 
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Figure SI: Experimental and computational results of the dynamic behavior of the 
system (associated to Figure 1). 

(A) Effect of the salicylate concentration in the growth rate of the cells (wild-type 
system, strain IE01). Absorbance values are averages of replicates. 

(B) Dynamic response of the system upon induction with different concentrations of 
salicylate (wild- type system). We represent the normalized fluorescence (YFP) with 
time. Error bars represent standard deviations. 

(C) Dynamic response of the wild-type system upon induction with different 
concentrations of salicylate in aerobic (dashed lines) and anaerobic (solid lines) 
conditions. We represent the absolute fluorescence (YFP) with time. Fluorescence 
values are averages of replicates. Anaerobic conditions were achieved by using 
mineral oil to cover the cultures. The no expression of the mar circuit with salicylate 
under anaerobic conditions may be explained, at least in part, by the fact that copper 
appears to shift intracellularly from Cu 2+ to Cu + in absence of oxygen (Rensing & 
Grass, 2003). 

(D) Computational simulations of the dynamic behavior with toy parameter values. 
Model-based sensitivity analysis of various parameters. (I) /? = 30 / 10 / 1. (II) p = 5 I 
2/1. (Ill) jto = 1 / 0.4 / 0.2 min" 1 . (IV) \i = 0.02 / 0.01 / 0.005 min" 1 . (V) k = 0.01 / 
0.02 / 0.1. (VI) 6= 0.5 / 0.1 / 0.01 min" 1 . For parameters not specified, we took the 
following parameter values: 6 = 0.5 min" 1 , a = 0.01 min" 1 , jto = 1 min" 1 , p = 5, k = 
0.01,^=30,z = c = 0, 6fe= 0.1 mM, v s = 2, anda=0.01. 

(E, F) Simulations to study the effect of Os and a on (E) response time, calculated as 
the time to reach the 50% of the steady state value (tso) upon induction with [Sal] = 5 

21 



mM, and (F) half maximal effective concentration of salicylate (Ks). We took the 
same parameter values as before. These plots show that when a approaches to 0 the 
response time increases because of the lack of active repressor in the cell. In the limit 
(a = 0), for [Sal] » &s, the situation will be like a gene without autogenous 
regulation. However, this is not what we observe experimentally, which leads us to 
consider a non-negligible value of a (e.g., a = 0.01). The plots also show that &s is the 
primary parameter responsible of changing Ks, although Ks can also be influenced by 
a in certain regime (e.g., 0.001 < a < 0.01). 



(G) Transfer function of copper regulation relating the fraction of oxidized MarR 
(MarR 0X ) with the relative amount of Cu 2+ . Experimental values (small circles, data 
taken from Hao et al, 2014) represent averages of replicates of the relative to the 



maximum oxidation. The line corresponds to 



1 r Cu 2+ 



-, fitted with K 0 



([MarR tot ]J 



1.35. At first order, we have a « 



K, 



[MarRj) 



v ° x [Cu 2+ ] 



. The intracellular copper levels 



upon induction with salicylate may be about 10 \iM (Rensing & Grass, 2003), and the 
concentration of MarR about 2 (Martin & Rosner, 2004). This gives a ~ 5%. 
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Figure S2: Experimental results and computational simulations of the dynamic 
response including rob and marB genes (associated to Figure 2). 

(A) Scheme of the natural circuit where MarR, MarA, MarB and Rob regulate the 
marRAB promoter. 

(B - D) Model predictions of the dynamic response (simulations to study the effect of 
Rob). (B) z = 100 simulates the wild-type system, and (C) z = 0 the Arob system. The 
arrows mark different response abilities at short times (blue curves). (D) Simulations 
(solid lines) together with gene expression experimental data (circles). We took the 
parameter values shown in Table S 1 . 

(E) Simulations to analyze the dynamics of the system for stable or unstable MarA 
(parameter S) and efficient or inefficient positive feedback (parameter p). MarR in 
blue, MarA in red. Black line corresponds to growth-rate-based exponential curve, 
We took the same parameter values as before. 

(F - I) Nullclines of the system (phase diagram). Solid black lines correspond to 
nullclines at 5 mM salicylate, and dashed black lines to nullclines without salicylate. 
(F, H) Wild-type system. (G, H) Arob system. (F, G) Theoretical trajectories (solid 
red curves for simulations starting from the origin; dashed red curves from the 
equilibrium point in absense of salicylate). (H, I) Representation of the normalized 
fluorescence (YFP, a measure of MarR) and the inferred promoter activity (a measure 
of MarA) with time in the phase space (red circles), together with the empirical 
trajectory calculated with Eqs. (S2) and (S3) (red line). The concentrations of MarA 
and MarR were normalized by the steady state values at 5 mM salicylate. 
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(J - L) Normalized dynamic response upon induction with 5 mM salicylate of (J) the 
wild- type system (strain IE01), (K) the Arob system (strain IE02), and (L) the AmarB 
system (strain DB01). Error bars correspond to experimental data, whereas the blue 

and black lines to different exponential models, [\-e~ Xt ^ and (l-e~^ respectively. 

The red line corresponds to promoter activity calculated with the best- fit exponential 
model. Fluorescence values are normalized by the maximum to compare curves. The 
boost time of the system (r) here was defined when promoter activity is maximal, Eq. 
(S5). We obtained the following parameter values: (J) A = 2.42 h" 1 , m = 1.66, ji = 0.42 
h" 1 , and r = 0.23 h; (K) A = 2.81 If 1 , m = 2.99, ^ = 0.43 if 1 , and r = 0.42 h; and (L) A 
= 2.28 If 1 , m = 2.92, ^ = 0.48 If 1 , and r = 0.54 h. 

(M) Dynamic response of the AmarB system upon induction with different 
concentrations of salicylate. Dynamics shown with solid lines correspond to the 
AmarB system, whereas dynamics shown with dashed lines correspond to the wild- 
type system. In both cases, we represent the normalized fluorescence (YFP) with 
time. Fluorescence values are averages of replicates. 

(N) Characterization by qRT-PCR of the intracellular RNA concentration of gene 
marA for different strains (wild-type, Arob, and AmarB), with (5 mM) and without 
salicylate. We report the ratio between the expression of gene marA and a constitutive 
housekeeping gene (16S ribosomal RNA, rRNA). Error bars correspond to standard 
deviations. 
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salicylate (mM) 



Figure S3. 
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Figure S3: Experimental results and computational simulations of the dynamic range 
and response time (associated to Figure 3). 

(A) Model prediction of the transfer function of the system. Relative promoter activity 
(in steady state) of the system as a function of salicylate. Solid black line represents 
model simulations, whereas circles (wild-type) and triangles {Arob) are experimental 
data. We took the parameter values shown in Table SI. In addition, the Hill-like 

TL min + ([Sal] / K Y s 

model H mar = — — — — — adjusted to all experimental data (from wild-type 

l + ([Sal]/K s ) ns 

and Arob systems; red line), gave U^ r = 0.10 (+/- 0.04), K s = 0.79 (+/- 0.10) mM, 
and ns = 1.36 (+/- 0.21). The plot shows comparatively theoretical and empirical 
values of input and output dynamic ranges. 

(B) Model predictions of the transfer functions of various control systems. We 
represent the relative promoter activity (in steady state) of the system as a function of 
the inducer. The blue curve corresponds to the wild-type mar circuit, whereas the red 
curve corresponds to a sensor gene negatively autoregulated and the yellow curve to a 
sensor gene constitutively expressed. For simplicity, here we took Vs = 1 and a = 0. 
The rest of parameter values are from Table S 1 . 

(C) The positive feedback strength (parameter k) and Rob concentration (parameter z) 
modulate the steady state of the system (simulations for 5 mM salicylate). Despite, the 
output dynamic range is almost constant. 

(D) Model predictions for different conditions (genotypes and salicylate levels). The 
predictions show Jt mar (5 mM)/;w(0 mM)| W T « ^w(5 mM)/;r mar (0 mM)| ra ^ « <x 2/3 = 
7.3, and Jt mar (5 mM)/;w(0 mM)\ marA « p = 10. In addition, they show 
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Jt m ar(WT)/jt mar (rob)\o mM « Jt mar (WT)/ jt ma r(rob)\ 5 mM « p « 2. Finally, they show 
JT mar (marA + )l ji mar (marA)\ marR « 1.4, and Jt mar {marA + )l jt mar {marA)\ marR> rob « p = 10. 
These predictions are in tune with previous experimental results (Martin & Rosner, 
1997). 

(E) Scheme of the coupled dual autogenous regulatory systems by CRP and MarR- 
MarA in the bacterial genome. The implementations are different each other. 

(F) Simulations to analyze the response time of two implementations of the dual 
autogenous control system. The blue curve corresponds to the wild-type mar circuit. 
The gray curves correspond to a hypothetical situation where the oxidized MarR acts 
as an activator. Solid gray curve for competitive binding between the repressor and 
activator, and dashed gray curve for independent binding. Parameter values from 
Table SI. 
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Figure S4. 
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Figure S4: Experimental results and computational simulations on the marRAB 
promoter integration of salicylate and cAMP signals (associated to Figure 4). 

(A, B) Dynamic response of the system upon induction with different concentrations 
of salicylate and cAMP. (A) Wild-type system (strain IE01). (B) Arob system (strain 
IE02). We represent the normalized fluorescence (YFP) with time. Fluorescence 
values are averages of replicates. 

(C, D) Normalized dynamic response of the wild-type system (C) and the Arob 
system (D) upon induction with 5 mM salicylate and 20 mM cAMP. Error bars 
correspond to experimental data (we repeated the experiment of induction for the 
wild-type system different times, but we always obtained high variability between 

l-e J 

and (\-e~^ respectively. The red line corresponds to promoter activity calculated 

with the best-fit exponential model. Fluorescence values are normalized by the 
maximum to compare curves. For the wild-type system, we got A = 1.60 h" 1 , m = 1.40 
and fji = 0.36 h" 1 , and for the Arob system, A = 2.50 h" 1 , m = 3.23 and fi = 0.43 h" 1 . The 
boost time of the system (r), here defined when promoter activity is maximal, Eq. 
(S5), was r = 0.26 h for the wild-type system, and r = 0.52 h for the Arob system. 

(E, F) Model-based predictions of the transfer function of the system as a function of 
the two input signals, salicylate and cAMP. (E) Wild-type system. (F) Arob system. 
For this plot, we used the parameter values shown in Table S 1 . 
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Figure S5: Theoretical analysis of the non-induced dynamic response (associated to 
Figures 5). 

(A - C) Phase diagram of the system in absence of salicylate showing an oscillatory 
pattern when MarA is cooperative (A). Here, we assumed high cooperativity (Hill 
coefficient of 4) for both MarA and MarR. The black lines correspond to the 
nullclines, with a different shape than for the wild-type system. The red curve 
represents a deterministic trajectory. The dynamics is represented with time in (B), 
where MarA is shown in gray and MarR in black. In (C) we show the resulting 
distributions of expression. For these plots, we took the following parameter values: 
6= 0.15 mm 1 , /a = 0.01 mm 1 , jto = 0.2 mm 1 , p = 10, 1, /3 = 5, andz = c = 0. 

(D - F) Phase diagram of the system in absence of salicylate showing stochastic 
pulses when MarA prevents MarR binding (D). To model this hypothetical 
competitive binding between MarA and MarR, we replaced yo in Eqs. (SI 1) by 
yo/(l+x). The black lines correspond to the nullclines, with a different shape than for 
the wild-type system. The red curve represents a stochastic trajectory. The dynamics 
is represented with time in (E), where MarA is shown in gray and MarR in black. In 

(F) we show the resulting distributions of expression. For these plots, we took the 
parameter values shown in Table SI, with z = 0. 

(G) Comparison of the distributions YFP and YFP/CFP for all single cells at all time 
points. The coefficient of variation (CV) is lower in case of YFP/CFP. 

(H) Scatter plot between YFP and CFP for a single cell trajectory showing a good 
correlation (dominance of extrinsic noise). 
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Figure S6: Experimental results of the dynamic response upon induction with 
salicylate at the single cell level (associated to Figure 6). 

(A) Box plot of the dynamic response (using YFP/CFP intensity) of multiple single 
cells for the wild-type system (strain IE01) and the Arob system (strain IE02). We 
present results for non-induced and induced (5 mM salicylate) cells. 

(B, C) Flow cytometry analysis (using YFP intensity) of the (B) wild- type and (C) 
Arob systems, for several conditions of salicylate. 

(D) Model simulations to study the stochastic response. Different trajectories 
(simulating different single cells) are shown in black. In green, we represent noise in 
gene expression (coefficient of variation). In blue, we show the dynamics in the 
deterministic regime (which approximates very well to the average). 

(E) Model predictions of noise in expression of MarA (r/ x , gray bars) and MarR (rj y9 
black bars) as a function of salicylate. 

(F) Experimental noise calculated with flow cytometry data, for the wild-type system, 
to validate the results reported in (E). 

(G) Map of the response time, which is determined by /) the time lag to react 
(empirical parameter m), and ii) the speed at which the expression increases 
(empirical parameter A). Open circles correspond the biological systems. The inset 
shows the experimentally determined dynamics. 

(H, I) Monitoring of fluorescence with time of four representative lineages showing 
different response times. (H) Wild-type system. (I) Arob system. 
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(J, K) Distributions of the empirical parameters A and m, obtained by fitting the 

\-e J to the single cell trajectories of 100 lineages. Black bars 

correspond to the wild-type system, and red ones to the Arob system. A and m are 
significantly higher in the Arob system (P < 0.05, [/-test). 

(L) Scatter plot of A and m [data from (J, K)] showing higher variability in case of the 
Arob system. 
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Figure S7. 
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Figure S7: Analysis of the dynamic response of a circuit with dual autogenous 
regulation. 

(A) Scheme of the dynamic response. 

(B) The response time depends with a decreasing trend on the negative feedback 
strength (understood as the binding affinity of the repressor to the promoter), being 
constant the maximal production rate of the operon (isoMPR). The inset shows model 
simulations. 

(C) The response time depends with a decreasing trend on the positive feedback 
strength (understood as the binding affinity of the activator to the promoter). The inset 
shows model simulations. 

(D) The response speed and response delay are positively correlated and they depend 
with an increasing trend on the degree of multimerization of the activator 
(cooperativity) and its ability to stimulate of RNA polymerase (fold-change). The 
inset shows model simulations. 
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Supplementary Tables 



Table SI: Nominal parameter values for the mathematical model. 



Parameter 


Value 


Note 


Reference 


5 


0.5 min" 1 


Protein half-life of ~1 min. 


Griffiths al, 2004 




0.007-0.01 
min" 1 




This work 


Jto 


0.1 min" 1 


Assuming transcription rate of ~1 
nM/min, translation rate of -0.05-0.1 
min" 1 , and mRNA degradation rate of 
-0.1 min" 1 


Levine et al, 2007 




30 


With RBS calculator 


Salis et al, 2009 


P 


10 


Experimental data show -6-fold 
increase in expression due to the 
direct effect of MarA. 


Martini al, 1996 
Chnhiz pt al 2012 


K 


0.02 


Having — 1UUU-1 jUU nJVl and 
^ = 5-10nM(as k = K r I K A ) 


Martm et al, 2008 
Seoane et al., 1995 


ds 


0.13 mM 


Adjusted to get the experimental 
transfer function. Estimated from K s 
~ 0.9 mM, knowing that 

0 s ^K s (2a) 1/vs 


Cohen al, 1993 


vs 


1.4 




This work 


a 


0.05 


Assuming [Cu 2+ ] / [MarR] ~ 5 upon 
induction with 5 mM salicylate 


Hao et al, 2014 


z 


100 


The number of Rob molecules per 
cell is -10,000, but they aggregate 
and are inactive. The fraction that is 
free is assumed -10% 


Griffiths al, 2009 


K c 


5mM 




Setty et al, 2003 


q 


0.3 


Knowing that q = r\ ex jt mar (2/ 


Swain et al, 2002 



38 



Table S2: E. coli operons with dual autogenous regulation 5 ' 6 from RegulonDB 



(Salgado et aL, 2013). 



Regulator 


Operon 


iNote 


Molecular 
implementation 


ChbR 


chbBCARFG 


Repressor that becomes 
activator in response to 
chitobiose. AraC/XylS family 
of transcriptional regulators. 


Multifunctional 
regulator 


CRP 


crp 


Repressor that becomes 
activator in response to 
cAMP. CRP-FNR 
superfamily of transcription 
factors. 


Multifunctional 
regulator 


LldR 


UdPRD 


Repressor that becomes 
activator in response to 
lactate. GntR family of 
transcriptional regulators. 


Multifunctional 
regulator 


MarR, 
MarA 


marRAB 


Control of multiple antibiotic 
resistance. AraC/XylS family 
of transcriptional regulators. 


Distributed 
regulation 



The operons with dual regulators (according to RegulonDB) ada-alkB, araC, dnaAN-recF, gadAXW, leuO, 
melR, murR, and phoPQ are not considered as systems with dual autogenous control in the sense of our study. 
6 AraC or MelR, also transcriptional regulators of the AraC/XylS family, repress their own expression, whereas 
they exert dual regulation (with arabinose or melibiose as inducers) on their targets. 
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